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1.  INTRODUCTION 


la  previous  papers  (Iglehart  and  Shadier  (1978a),  (1978b),  (1979) 
and  Shadlar  (1979)),  .-a  provided  methods  for  obtaining  from  a  single 
simulation  run  point  eseiaaces  and  confidence  inCervais  for  general 
characteristics  of  "passage  times"  in  cercain  closed  networks  of  queues. 
Informally,  a  passage  ciae  is  the  time  for  a  job  to  traverse  a  portion  of 
the  netvork.  Such  quantities  are  important  in  computer  and  communication 
system  models,  and  In  this  context,  quantities  other  chan  mean  values  are 
of  Interest,  the  basis  for  these  estimation  methods  is  the  regenerative 
method  for  simulation  analysis  (Crane  and  Iglehart  (1975)).  For  an 
introduction  to  and  a  detailed  review  of  the  regenerative  method,  see 
Crane  and  Lemolne  (1977)  and  Iglehart  il^75a). 

We  consider  here  the  calculation  of  theoretical  values  for  variance 
constants  entering  into  the  central  limit  theorems  used  to  obtain 
confidence  intervals  from  passage  time  simulations.  Using  results  of 
Hordijk,  Iglehart,  and  Schassberger  (19*6)  for  the  calculation  of  moments 
in  discrete  time  and  continuous  time  Markov  chains,  we  calculate  variance 
constants  pertinent  to  mean  passage  times.  We  do  this  first  for  the 
broadly  applicable  "marked  Job  method"  for  passage  time  simulation  vhlch 
is  based  on  the  cracking  of  a  distinguished  Job.  and  thea  for  the 
"decomposition  method"  in  which  observed  passage  times  for  all  of  the  jobs 
enter  into  the  construction  of  point  and  interval  estimates.  The 
decomposition  method  provides  point  and  interval  estimates  for  a  restricted 
but  practically  important  class  of  passage  times,  namely,  those 
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corresponding  Co  the  passage  through  a  subnetwork  of  the  given  network  of 
queues . 

The  results  of  this  peper  provide  a  firs  bAsls  for  compering  the 
statisticsl  efficiency  of  the  cvo  methods  vh«a  boch  Apply.  The  oaIcuIacIoqa 
Also  permit  aq  assessment  of  the  AfflCAcy  of  the  narked  job  method  for 
simulation  of  "response  cia«s"  (ccoplsct  circuits  ia  the  network) ;  the 
marked  Job  method  is  Apparsaely  the  only  available  means  of  obcsiniag 
coaiidsacs  intervals  for  r«spoas«  class  from  a  siagls  siaulACioa  run. 


2.  ?*£LlMIlua.:SS 

Vs  ooasldsr  olossd  networks  of  queues  with  a  fiaics  number  of  Jobs 
(customers) ,  M.  la  each  as ever  it  there  art  a  fiaics  nuabsr  of  service 
esaesrs,  s,  aad  a  fiaics  nuabsr  of  Job  olassss ,  c.  Ac  aach  epoch  of  clas 
each  Job  is  ia  exactly  oas  Job  class,  bu;  Jobs  may  chaags  class  is  chsy 
traverse  chs  asevork.  L’poa  ccoplecicn  of  service  at  ceacer  i,  a  Job  of 
class  J  joss  to  ceacer  k  md  changes  co  class  1  with  probability  ?  4  ,  Vs 

e»  —  •  Na, 

assuae  chac  ,  .  :lii,kss,  l£J,iSc;  is  a  given  irreducible  Markov 

L  j  * 

aacrix. 


AC  eaca  service  cancer  Jobs  queue  and  receive  service  according  co  a 
fixed  priority  scheme  aacng  olassss,  which  scheme  can  vary  iron  center  co 
center ■  Each  center  operates  as  a  single  server,  processing  joes  of  a 
fixed  class  according  to  a  fixed  service  discipline.  All  service  tines 
in  the  network  are  mutually  independent,  and  at  each  cer.tar  have  a 
distribution  with  a  Cox-phase  (exponential  stage)  representation 


i 
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(cf.  Cox  (1935).  Galanba  and  hunts  (.1976))  with  paramacars  which  nay  dapand 
on  tha  sarvica  cantar,  class  of  job  balng  sarvlcad,  and  cha  "stats"  of 
cha  antira  systaa.  (Wa  axcluda  taro  sarvlca  ciaas  occurring  vich  posltiva 
probability.)  A  Job  in  sarvica  aay  or  aay  noc  ba  praaepcad  (according  co 
a  fixad  procadura  for  aach  cancar)  if  anochar  job  of  highar  priority  Joins 
tha  quaua  at  cha  cantar.  Wa  rastrict  cha  prasanc  discussion  co  nacvorks 
in  which  all  sarvica  cimas  ara  axponandally  distrlbutad,  and  daal  with 
distrlbucions  having  a  Cox-phasa  raprasancation  in  cha  usual  way  by  cha 
nachod  of  stagas. 

As  in  Iglaharc  and  Shadlar  (1973a),  wa  vlaw  cha  S  Jobs  as  baing 
cooplacaly  ordarad  in  a  linaar  scack,  and  darina  cha  vaccor  Z(c)  according 
to 

2(c)  -  (c[1)  (t) . C,(1)(t),S.(t);...; 

\  Jk(l)  J1  1 

C.(,)  (c) . C4(#)(t),S  (t)\  .  (2.1) 

Jk(s)  n  *  / 

Tha  linaar  scack  corrasponds  co  cha  crdar  of  ooaponancs  in  cha  vaccor  Z^c) 
afcar  Ignoring  any  cars  ooaponancs.  Within  a  class  ac  a  cancar.  Jobs 
waiting  appaar  in  cha  linaar  scack  in  cha  crdar  of  chair  arrival  ac  cha 
cancar,  cha  lacasc  co  arriva  balng  closasc  co  cha  cop  of  cha  scack. 

Lattlng  N\c)  danoca  cha  position  froa  cha  cop  of  cha  aarkad  Job  in 
this  linaar  stack,  for  ciO  cha  scats  vaccor  of  cha  nacvork  is 

XU)  -  (2(0, NU))  .  (2.2) 
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As  bafara,  wa  apaciiy  tha  paasaga  tlaa  for  tha  sarkad  Job  by  four  subsats  A,, 

A,,  3^,  and  3,  of  £,  tha  aut«  spaca  of  tha  procaaa  5«{X(c> : :*0> .  Tha  aaca 

A,  and  A,  [raap.  3.  and  dacaraiaa  whan  to  start  [raap.  atop’  tha  clock 

aaaaurlag  a  particular  paasaga  tiaa  for  tha  sarkad  Job.  Danociag  tha  juap 

ciaaa  of  X  b*  (t  :uiO},  for  k.ail  wa  raqulra  that  tha  aaca  A,,  A, ,  3,  and 

3,  aaclafy: 

* 


and 


if  XCt^JcA,^,  X(tq)«A:,  and  XCr^aA,  , 

than  X(r  ,  )«3,  and  Xyt  '*3,  for  scan  >ask. 

3*  i  rx^in  * 

if  X<ta'«3,,  and  Xv^.S,  . 


n-.-tc 


S*K' 


than  Xvt  ,  )«A,  and  Xyt _ 'aA,  for  ac®a  Tsa'k. 

a*  i  O'*®  * 


Thaaa  ccndiciona  anaura  that  tha  atart  and  tarainacicn  ct_sas  for  :ha 
apaciiiad  paasaga  tlaa  strictly  altarnata.  Also  la  tarns  of  thaaa  Juap 


tlsaa , 

va  dafiaa 

two  aaquaacaa  of 

Tha  at 

art  [rasp. 

tamination]  tlaa 

Job  la 

ianctad  by 

3t  ,  [rasp.  T4’. 

-  “*  4 

sarkad 

Job  bagisa 

at  t»0,  wa  hava 

*0 

•  0 

3, 

4 

•  iaf  { niT,  ::Cv*a' 

a* 

*  4 

•  iaf  a>S,_. ;Xyt 
*  ■* 

n-. 


Than  tha  jth  oassaga  tisa  for  tha  sarkad  job  la  ?, •?, -i ,  .  .  Jii.  N'ota 

^  ^  a 

that  thaaa  iafiaiticr.s  ara  aoaawsat  sera  raatrictiva  than  thcia  ir.  [10] 
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W*  let  Xn  denote  the  ium  of  th«  process  X  whsn  the  (n-t-l)sc  passage 
ciae  of  chs  marked  job  begins:  X  -X(S J ,  niO.  The  process  {(X  ,P  ^.):niO} 

u  a  Q  n+L 

is  a  regenerative  process  in  discrete  tine,  and  the  regenerative  property 
guarantees  (Miller  (1972))  that  as 

<Xa’W  <X’P>  • 

where  — >  denotes  convergence  in  distribution.  The  random  variable  ?  is 
the  limiting  passage  time  for  the  marked  job,  and  the  argument  in  the 
Appendix  of  [9]  shows  that  the  sequence  of  passage  times  for  any  other 
job  also  converges  in  distribution  to  the  same  random  variable.  The  goal 
of  the  simulation  is  estimation  of 

r (f )  -  E(f (P) }  ,  (2.3) 

where  f  is  a  real-valued  measurable  function  with  domain  R+“[0,®).  We 
assume  that  ?{P£D(f)}«0,  where  0(f)  is  the  3et  of  discontinuities  of  the 
function  f. 

3.  THEORETICAL  VALUES  FOR  FINITE  STATE  MARKOV  CHAINS 

Ut  {X^:  k20}  be  an  irreducible  Markov  chain  with  f'  <te  state 
space  E*{0,1 . N}  and  one-step  transition  matrix 

P  -  (p^si.jaE}  . 

For  ail,  denotes  the  n-step  transition  probability  from  state  i  to 
state  j ,  and 
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For  a  fixad  suu  i«E,  ?^{*}  danocas  cha  condicionai  probability 
associacad  vlch  scareiag  cha  chain  in  aeaca  1,  and  E^*}  danocas  cha 
corraaponding  condicionai  axpaccacion.  For  j«E  aad  nil,  va  lac  Sa(J) 
daaoca  cha  nch  ancraaca  ciaa  of  {X,  :'*20)  :o  scaca  j,  a.g., 

3,  (J )  ■  aiai kii :  > 
aad  lac  and  3aU)*SaCi>*'Sa«.1C2)  •  a>l. 


Wa  caasidar  vaccars  such  a a  (v(Q) ,v(l) , . . . ,v(M) )  ca  ba  column  vaccars, 
and  '/lav  raai-vaiuad  functions ,  auch  as  f  and  g,  having  domain  I  in  chis 
way.  - nlaaa  apaciiiad  acharvisa,  cha  symbol  !(•;  danocas  cha  vaocar 


V  —  qV  *  4  >  E^l  ******  >  '  *  • '  ) 


la  addition  (far  vaccars  u  and  v)  cha  symbol  u*v  danocas  cha  vaccar 


(u(O)vi^O)  ,a(l)vvl) ....  ,u(N)v(.20) 


'or  a  nacrix  Am  . ,a  ),  va  lac 

'1  a.  a 


u*A  •  A»u  •  ([u»a-,uaa . u»a  ) 

J  •  s 


Finally,  for  a  nacrix  3*vb..,b . ,b  ),  va  la: 

J  .  a 


A*E  •  „a,»3,,a,  »b . a  »b  ) 

o  j  .  *.  a  a 


For  cha  discraca  ciaa  Markov  cnaia  ;  : kiO ; ,  va  oonsidar  hara  only 
oyolaa  of  cha  raganarici"a  procaas  foraad  by  cha  succassi *a  ancranoaa  :: 


scaca  3.  and  hancafarch  aupprasa  cha  0  in  cha  notation  c  vO) ,  n  vO) .  see. 

U  a 


Far  < I  aad  n«0,l 


lac 
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o?“]  *  W*-  V;!  • 

•ad  ««c 

0P  "  * 

We  obtain  ^P^P  fr°«  P  by  setting  ch«  0-column  of  P  equal  to  0.  It  is 
easy  to  s««  that  QPa  is  cha  matrix  product  of  n  copias  of  QP,  and  that 
for  all  ail.  0?J0“0* 


For  any  raal-valuad  function  f  with  domain  E,  cha  state  space  of 
{X^:kiO},  we  define 


**>  ■  £ 
•  I _ A 


Y  (f)  -  2.  fcx,  )  . 
k»0 


Theorem  (3.1)  of  Hordijk,  Iglehart,  and  Schassbergar  (19’6)  shows  that 
for  an  irreducible,  finite  state  discrete  time  Markov  chain  with  transition 
matrix  ?, 


E{Y1(f)>  -  (I-0?) 


(3.1) 


E(Yi(f)Y1(g)}  -  a-0?)"^h  , 


(3.2) 


where  h-f sEtY^  (g>  }<*g®StY,  (f)>-f»g. 


Now  we  consider  continuous  time  Markov  chains  JC"{X(t) : CiO;  with  finite 

state  space  E-{0,1 . N},  transition  matrix  ?(c)*{p^ (c) : i,J«E)  and  matrix 

of  infinitesimal  parameters  Q-iq^ : i,J<S} .  The  exponentially  distributed 


a 


holding  time  is  any  scat*  icE  ha*  rat*  parameter  q^«-q^  .  For  all  lc£,  w* 
assume  that  0<q1<»,  1.*.,  that  all  stacts  art  stable  and  nonabsorbing , 
and  la  addition  that 


Z  < 


i-0 


ij 


0  . 


This  last  assumption  guarantees  that,  scartiag  from  any  state  i«E,  the 
Markov  chain  $  sake*  a  transition  to  a  next  state  j«E.  V*  now  tors  the 
Jump  matrix  E»{r^}  of  th*  chain,  defining  the  elements  r..^  according  to 

m 

a 

0  if  j-i 

Ve  assume  that  the  Jump  matrix  3.  is  Irreducible  vand  therefore  positive 
recurrent).  This  is  equivalent  to  the  continuous  time  Markov  chain  X  being 
irreducible.  Far  J«E  and  nil,  we  let  3  (J)  denote  the  nth  entrance  time 
a:  X  to  state  J ,  i.e. , 


3,  (J)  *  iaits>0:X«.s-)»J  ,  X^s)' 


As  la  the  case  of  discrete  time  Markov  chains,  ve  restrict  attention 
to  regenerative  cycles  formed  by  the  successive  entrances  to  state  3,  and 
suppress  the  3  in  our  notation.  For  taO,  ve  let 


and.  for  niC.  construe 

constructed  ,?a  from  ? 

0 


the  aatrix  ,3.“  from  3  in  the  same 
in  the  discrete  time  rase. 


as  ve 
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vita  scare  space  l.  3enocing  by  L(c)  the  last  state  visited  by  the  Markov 
chain  X  before  jumping  to  X(c),  the  process  V«{V(t):taO;  where 

VCt)  -  (L(t).XCt))  , 

is  the  fundamental  stochastic  process  of  the  passage  tine  simulation. 


Ihe  process  7  has  a  state  space,  F,  consisting  of  all  pairs  of  states 
i.j«S  for  which  a  transition  in  X  from  state  i  to  state  j  can  occur 
vi th  positive  probability.  Since  X  is  an  irreducible,  positive  recurrent 
Markov  chain,  30  is  V.  Ve  define  subsets  3  and  T  of  the  state  space  F 
according  to 


and 


3  *  {(k,3)sF:k£A1 ,  aeA,; 


I  *  { (k.a) £ J : ks3, ,  ae3,}  . 


The  entrances  of  7  to  the  fixed  subset  3  [rasp.  Z]  correspond  to  the  starts 
"rasp,  terminations’  of  passage  tines  for  the  narked  job. 


As  in  Iglahart  and  Shadier  (1973a' ,  we  select  a  (fixed 


j  state  or  a. 


designated  state  3,  and  assume  that  7(3)*0.  To  estimate  r(f),  the  narked 


jco  nethod  prescribes  that  we  carry  out  the  sinulat 
defined  by  the  successive  returns  to  state  3;  with! 
the  number  or  passage  tines  of  the  narked  job  and  n 


;icn  of  7  in  3-cycles 
,n  each  cycle  we  record 
easure  each  of  these 


passage  tines 


LI 


We  let  {V  :n20}  denoce  Che  embedded  jump  chain  associated  with  the 
continuous  time  process  V.  The  random  times  {y  ; nil }  denote  the  lengths 

tl 

in  discrete  time  units  of  the  successive  0-cycles  (successive  returns  to 

the  fixed  state  0)  for  {V  :n*0},  and  we  define  (5^*0  and  5  ■y.+.-.+y  ,  mil. 

ti  u  ml  m 

Then  the  number  of  passage  times  for  the  marked  job  in  the  first  0-cycle 
of  V  is 

sr1 

Ml  "  i?o  l(Ysl  ' 


(For  a  sec  A,  1  (x)*l  if  xeA  and  0  if  x^A.  Here  we  suppress  the 
argument  u.)  The  sum  of  the  values  of  the  function  f  for  the  passage  times 
of  the  marked  job  in  that  cycle  is 


Yx(f) 


2J  £(?,)  . 

j-1  J 


We  denote  the  analogous  quantities  in  the  kch  0-cycle  by  and  Y^Cf") 


The  key  results  leading  to  point  estimates  and  confidence  intervals 
for  r(f)  are  that  the  pairs  of  random  variables 

{(Y^CfJ.iy  :kil}  (4.1) 

are  Independent  and  identically  distributed,  and  that 

r(f)  -  E0{Yl(f)>/E0t'M1>  ,  (4.2) 


provided  chat  the  quantity  E{ j f (P) [ }<«. 


r 


L 
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Given  Equations  (4.1)  and  (4.2),  eh*  regenerative  sec hod  provides  from 
a  fixed  number  a  of  0-cycies  the  so-called  classical  point  estimate  (cf . 
Iglehart  (1975b)) 


r  (f)  -  Y  (f )  /M  , 

Cl  Q  Q 


and  the  associated  confidence  interval  for  r(f)  follows  from  the  central 
limit  theorem 


LI/2{rQ(f)-r(f)} 

a(t)/E0(ML; 


->  51(0,1) 


(4.3) 


Here  51(0,1)  is  a  standardized  (naan  0,  variance  1)  normal  random  variable 
and  z"(f)  is  the  variance  of  Y, (f )-r (f )M7 . 

"or  calculation  of  theoretical  values,  we  restrict  attention  to  the 
nean  passage  time;  thus,  the  function  :  in  the  definition  of  r(i)  is  the 
identity  function.  Using  the  results  of  Section  2.1,  ve  shew  hew  to  compute 
the  value  of  the  mean  passage  time  r  and  the  corresponding  variance 
constant  r2  appearing  in  the  central  limit  theorem  of  Equation  (4.3).  These 

calculations  rest  on  the-  definition  of  two  particular  functions  ^denoted 

*  * 

f  and  g  )  having  domain  F  and  taking  values  in  the  set  -.3,1;. 


Ve  define  the  function  f  to  be  tne  indicator  function,  1. ,  of  the 
sec  S;  i.e.,  for  (z ,n. z ' ,n' )<?, 


-  .  -  «. > 

•  f  / 
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Proposition  4.1  follows  directly  from  Equation*  <3.1)  end  (3.2) 


PROPOSITION  4.1.  Lsc  f  bs  ths  function  defined  by  Equation  (4.4),  «nd  R 
chs  trsnsltion  mstrlx  of  ths  dlscrsts  tins  Markov  chsln  {V^ikkO}.  Thsn 


rl'1  .  ) 

t(vf •  '\&  ‘  * 


vhsrs  i,  is  chs  elms  of  chs  first  rsturn  to  ths  ststs  0  in  (V,  :kiO)  and 
1  k 

•  *  *  ,  *  * 
h  -:f  •EiTl(f  )}-f  *f  . 


Vs  uss  this  rssult  and  chs  definition  of  M,  to  obtain  ths  quantities 


and  Ej(M‘}  according  to 


VMi}  •  vvfV 


E,(M?}  -  E.wT.vf  >)' 

0  1  0  A 


For  an  clamant  (a.n.a1 ,n' )«F,  ths  value  of  ths  function  g  is  1  if  a 
passage  else  for  chs  marked  Job  starts  or  is  underway  when  the 
configuration  of  chs  Job  tcack  is  i'  and  chs  marked  Job  is  in  position  n' . 
and  is  0  otherwise.  Formally,  1st  2  be  chs  scace  space  of  the  process 
£av2(e)t tiO},  and  denote  by  C  ths  set  of  veentsr .class'  pairs  ;n  the 
network.  Vs  define  a  function  h  taking  values  in  0  and  having 


domain 


n  1,2, . . .  ,50  as  follows.  For  2*2  *nd  Isnsh’,  ths  value  of  hv.:,n 


J 


u 


1*  (i.J)  vhan  cha  Job  la  position  a  la  eh«  job  stack  t  la  of  cl*aa  j  at 
cancar  1.  Nov  considar  cha  ambaddad  Jump  chain  {V^:kiO}  aaaoclacad  with 
cha  coaclnoua  elaa  Markov  chain  V.  For  scacas  v',v"cF,  cha  scata  apaca  of 
{V^skiO},  va  vriea  v’-vv*'  vhaa  v"  la  accaaalbla  from  v',  l.a.,  vhan  for  aoma 
nil ,  cha  probability  *  carting  from  v'  of  an tar lag  v"  on  tha  nth  atap  la 
poalciva.  Similarly,  for  any  aubaat  L  of  F  va  vrlct  v’  /\>  v"  vhan  v"  la 
accaaalbla  from  v'  undar  tha  taboo  L.  This  naans  (cf .  Chung  UW?),  pp .  »5, 
*&)  chat  for  aoma  nil,  chart  la  a  poalciva  probability,  starting  from 
scata  v' ,  cf  ancarlag  v"  oa  tha  nth  stap  undar  tha  rascriction  that  nona 
of  cha  scacas  la  L  Is  aacarad  la  bacvaan  vaxduaiva  of  both  ands). 


Canoeing  tha  sac  of  (cancar, class)  pairs  in  tha  naevork  by  C,  va  dafina 
a  subaac  G  of  C  according  to 

G  •  {(i,J)«C:  for  soma  (8, a, s' ,a’ )«S,  h(* ' ,n* )•(!, j ) }  u 

;(i,J)<C:  for  soma  (s.a, s' ,n' )<F-(SuT) ,  v’*S  and  v"<T  , 

T  “ 

.  -  .  <  -is 


la  indarsay] 


g  v8 , a, 8 ' .a’) 


<.«v"  and  hi, a 

’  ,  a  ’ )  ■  ( 1 ,  j  )  ;  . 

'.class)  pair 

is  in  tha  sat  G,  [ rasp . 

4k 

ob  to  ba  of 

class  j  ac  cancar  1  vhan 

'  A,  ,  ,  3.  , 

L  a»  Jk 

and  3.  starts  [rasp. 

k 

t^.t  function 

i  » 

la  •  *  a'U 

A  » i  \  *  *  U  * 

L 


IS 


Than  v«  have 

PROPOSITION  4.2.  Lac  g*  ba  cha  function  defined  by  Equation  (4.7),  and  R 
ba  the  Jump  matrix  and  q  tha  vactor  of  rata  paranatara  for  bolding  times 
in  tha  continuou*  tima  Markov  chain  V.  Than 

E{Yl(g*)}  -  e|J^  1  g*(V(a))dsj  -  (I-QIl)’l(g**q“l)  , 

and 


EtCY^g*))2}  -  (l“0R)”1(h**q_^)  , 

whara  8,  la  tha  tiaa  of  tha  flrat  raturn  to  tha  stata  0  in  V,  and 
l 

h*-2f**E{YL(g*)}. 


Proposition  4.2  follows  diraccly  from  Equations  (3.3)  and  (3.**) .  Ve 
usa  this  raault  togathar  with  tha  observation  that 


g*(V (s) )da 


V».3) 


to  obtain  tha  quantities 


and 


.  o\ 


C- . 10) 


'J«la*  tha  r»da  foraula,  Equation*  (4.5)  and  (i.9)  yiaid  tha  quantity  r. 
To  obtain  an  axprasaion  for  tha  varianca  cooataat  «rfc  appaaring  in  tha 
•ioocrai  Liait  thaoraa  (Equation  (4.3))  for  cha  aarkad  Job  aathod,  wa  raquira 
ona  additional  raault. 
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Whaa  comparing  tha  statistical  afficiaacy  of  ch«  markad  job  and 
dacomposition  mat hods,  It  is  convaniant  to  h«va  a  cancral  Halt  thaoram 
comp Arab la  to  Equation  (4.3)  but  la  tarma  of  simulation  tint,  c,  rachar 
chan  nuab«r  of  cyclas,  a.  Lac  a(c)  bA  chA  numbar  of  passaga  timas 
coaplatad  by  tlma  t,  I.A.,  la  ChA  latarval  (0,t].  If  wa  dAaoCA  by  alt)  thA 
nuabAr  of  0-cyclas  complatad  by  tima  t,  thAa  from  r anawal  thAory,  as  t-**», 


£i£l 

c 


with  probAblllcv  oaA,  vhara  E-(a,}  la  chA  axpActAd  langth  of  a  0-cyclA  In 
V.  This  Imp  Has  that  foe  Largs  t,  thA  nuabAr  of  0-cyclas  oomplatAd  by 
Clam  t  la  approxiaataly  e/E^ia^}.  Combining  this  raault  with  Equation  (^.3), 
lc  follows  chAC  aa  c-*», 


1  /  *? 

.,/«(c)  \ 

c1. 

|U(t)>  ^  iK?L)J  -rCf). 

a 

HlO.i)  . 


This  raault  is  indapandant  of  tha  initial  scaca  V(0).  Sines  tha  numarator 
la  this  cantral  limit  thaoram  la  indapandant  of  tha  stats  0  salactsd  to 
form  oyclaa.  so  is  tha  danominacor .  Thus  for  tha  maan  passaga  elms. 


A 


«Ai|W,,W 


is  chA  appropriaca  maasura  of  stacistioal  affioisney  for  tha  aarkad  job 
aathod  and  is  indapandanc  of  tha  scats  0«S  aaiactad  to  fora  eye  las.  Nets 
chat  wa  obtain  tha  quantity  according  to 


YLa>} 


* 


W  • 
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where  1  is  the  £  function  identically  equal  to  one  end 


E{Y1(1)} 


sffl 


l(V(s))ds 


(I-QR)'LCl»q‘1) 


5.  VARIANCE  CONSTANTS  FOR  THE  DECOMPOSITION  METHOD 

W«  now  turn  to  the  decomposition  method.  As  la  Shedler  C1979)  we 
label  the  jobs  from  1  to  N,  and  for  1*1, 2,..., N,  denote  by  N^(t)  the 
position  of  job  1  In  the  linear  job  stack  at  time  t.  Then,  la  terms  of  the 
vector  Z(t)  corresponding  to  the  job  stack,  we  set 

X°(t)  -  (2(t),.'TLCt),N2Ct) . N*\t))  . 

0-0 

The  process  X  *\X  ^t):tiO;  is  an  irreducible,  positive  recurrent  continuous 
else  Markov  chain  with  state  space  . 


We  let  L°(t)  denote  the  last  state  visited  by  the  Markov  chain  3?  before 
jumping  to  X^(t),  and  for  t20  define 

V°(t)  -  (L0(c),X°(t))  . 

The  process  7^*{7(t) : tiO;  is  the  fundamental  stochastic  process  of  the 
passage  time  simulation. 
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T  *  { (t| •  •  •  . S  t  *  *  * i tl^|)  SF  •  f or  IOBO  £ | 

(z.n^)*^,  end  (s’,a£)cB2}  . 

The  successive  entrances  of  V3  co  eh*  fixed  subs«c  S°  of  F°  correspond  to 
eh*  scares  of  passage  class  irrespective  of  job  ld*ncley,  and  eh*  entrances 
of  V13  to  eh*  subs*e  T13  correspond  eo  eh*  terminations. 

Th«  decomposition  a*chod  applies  to  passage  tia«s  for  which  th*  sees  S 

and  T  (which  d«fin*  th*  scares  and  e«rainacions  of  passag*  tines)  ar* 

disjoint.  D«nocing  by  {P3:nil}  eh*  sequence  of  pass ag*  tines  (irrespective 

of  Job  identity) ,  enumerated  in  order  of  passage  time  start,  by  the 

argument  in  eh*  Appendix  of  [9],  ?3->P° ,  and  the  goal  of  the  simulation 

a  a 

is  estimation  of 

r°(f)  -  E(f(P°)>  . 

To  estimate  r3(f) ,  we  carry  out  the  simulation  of  the  process  V*3  in 
random  blocks  defined  by  th*  successive  entrances  of  the  process  to  the 
fixed  set  of  scares  U3 .  Entrances  of  V3  to  th*  set  'J3  correspond  co  the 
terminations  of  passage  times  (irrespective  of  Job  identity)  which  occur 
when  no  other  passage  times  are  underway,  and  which  leave  a  fixed 
configuration  of  th*  Job  stack.  Formally,  recall  that  C  is  the  set  of 
(center .class)  pairs  in  the  network  and  define  a  subset  H  of  C  according 

H  ■  ((i,J)«C:  for  some  (z.n.z’  ,n’  )«T-S,  ht.z '  ,n'  )«(i,j ) ;  u 
{(i.J)«C:  for  some  (z.a.z' ,n’)<F-($uT) ,  v'<T  and  v"«S  , 
v'  4s  (z.n.z’.n'),  (z.n.z'.n')  &  v"  and  h(z ' ,n’ )«(i,J ) }  . 


to 
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Thus,  the  sec  where  *  (center .class)  pelr  Is  In  the  sec  [resp. 

Hjl  if  it  Is  possible  for  the  marked  job  to  be  of  dess  j  at  center  i  when 
e  passage  tine  for  the  narked  job  specified  by  the  secs  A^,  A^,  B^,  and  3^ 
terminates  [resp.  is  not  underway]. 


Now  define  a  subset  0°  of  D,  the  state  space  of  Z»{Z(c) : ckO} ,  according 
CO 

•  {scO*. h(c.,n)<ff"fbr  ISqSN  and ^ for  some  n„.  h(r,a)«E^}  • 

Elements  of  the  sec  3^  correspond  to  configurations  of  the  job  stack  upon 
termination  of  a  passage  time  with  no  ocher  passage  times  underway.  The 
sec  3^  is  nonempty  because  SfiT*»,  and  thus  the  sec  3.,  is  nonempty. 
Therefore,  we  can  select  an  element  2°  of  33,  and  (in  terms  at  this  fixed 
2J)  define  the  set  according  to 


,0  ,  0, 

:  2  *2  ; 


t 


where  T^  is  the  subset  of  corresponding  to  the  terminations  of  passage 
times. 


0 

.-or  Veil,  y.  denotes  the  length  in  discrete  time  units  of  the  kth  block 
0  .  ^ 

(returns  to  the  sec  u  )  of  the  process  -.7  :aiQ;;  the  latter  is  the  embedded 

n 

jump  chain  associated  with  7®,  and  we  set  o?«0  and  c^«y^.  .  .^vV ,  mil.  Ve 

u  a  a  a 

0  0  0 

assume  that  V  (0)<1  and  denote  by  :<  the  number  of  passage  times  in  the 

a 

0  0 

meft  block  of  the  process  7  .  Also,  we  let  ?  (f)  be  the  sum  of  the 

a 

0  3 

quantities  i(?  )  over  the  passage  times  in  the  nth  block  of  7  .  The  key 

results  of  [1*J  leading  to  point  ascinaces  and  confidence  intervals  for 

the  quantity  r  (£}  are  that  the  pairs  of  random  variables 
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UY?(f).K°):ail} 


(5.1) 


are  independent  and  identically  distributed,  and  chac 

r°Cf)  -  *  o{Tl<f)}/E  0{*1}  .  (5.2) 

provided  chat  tha  quantity  £( f  (P°)  Tha  symbol  E  .{•}  is  an  abuse  of 

Uu 

our  pravlous  notation.  It  connotes  conditional  expectation  associated 
with  starting  tha  Harkov  chain  V°  in  one  of  tha  states  in  the  sec  U°.  The 
definition  of  the  sac  U°  implies  chac  the  conditional  expectation  is 
independent  of  the  particular  starting  scace  in  l'3. 


Given  Equations  (5.1)  and  (5.2),  from  a  fixed  number  of  blocks  of  V  , 
the  decomposition  method  provides  the  point  estimate 

r°(f)  ■  . 


and  tha  associated  confidence  interval  for  r  t>£)  follows  from  the  central 
limit  theorem 


n1/2{*jj(f)-r°(£)} 


N(0,1) 


The  quantity  (ju(£))*  is  the  variance  of  Y*?(f )-r^(f )K?. 

1  x 


(5.3) 


Taking  :  to  be  the  identity  function,  we  restrict  attention  to  the 
aaan  value  r ^  and  consider  the  associated  theoretical  values.  3v  the 


argument  which  leads  to  Equation  v*.12)«  an  appropriate  measure  of  the 
statistical  efficiency  of  tha  simulation  is  the  quantity 


I 


cs.i) 


where  is  the  Length  of  s  block  in  the  continuous  tine  process  v®. 


The  individual  quantities  required  to  compute  this  measure  of 
efficiency  sre  defined  in  terse  of  the  successive  returns  of  the  process 
to  a  fixed  sec  of  scsces  rsther  then  to  s  single  state.  Moreover, 
the  successive  entrances  of  V*3  to  ’O^  are  not  regeneration  points  for  V^. 
Accordingly,  it  is  noc  possible  to  apply  the  results  of  Section  1  directly, 
as  ve  did  for  the  narked  Job  nechod.  %'e  can,  however,  select  a  fixed 
state  ^designated  state  uu)  from  the  sec  1’^  and  compute  the  quantity 
corresponding  to  equation  (3.i)  for  the  resulting  u^3 -cycles .  v^oce  that 


te  successive  entrances 


of  the  process  to  the  fixed  st 


ate  u  are 


regeneration  points  for  the  process  7  .}  The  expression  in  equation  v.5.*) 
computed  for  u^ -cycles  is 

o  cN  _  o.,i  ;.o  .  ,-.o. 

*  u  )  •  \ c  Oil, :  7  i  Oi.v,  •  . 

'4  *  J  u  * 

2  j  0 

where  the  constant  ^analogous  to  ;  '  is  defined  for  u  -cycles .  This 

J  0  0  ^ 

quantity  e  vu  )  is  equal  to  e  .  To  see  this,  for  tiO  let  n\t'  :e  tne 

number  of  passage  tines  v irrespective  of  job  identity'  completed  in  the 

interval  ,3,t].  In  tens  of  tinulation  tine.  t.  ve  have  the  central  line: 

theorem 


'  d  a  •  *  • 

*  « 


cl 

.  *c  *  •  \  mb  K  >  M 

T*  *  i 


anc .  when  :  is  the  identity  Junction,  the  variance  constant  m  the 
denominator  is  the  quantity  e' .  There  is  a  similar  central  linit  tnercen 
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la  tars*  of  u^-cycles;  tha  numerator  la  cha  saaa  and  tha  variance  constant 
la  tha  danoainator  la  e^u^) .  Slaca  tha  auaaracors  in  chaaa  tvo  central 

limit  chaoreme  are  cha  same,.  a  a  ara  cha  Halting  random  variables  <N(G,1)), 
e^  auat  equal  a°(u'"'/.  For  a  similar  argumenc,  saa  Propoaltiona  3.1  and 
3.6  of  Crana  and  Iglahart  (1973). 

Naxt  ua  observe  chat  tha  numbar  of  passaga  t  la  a  a  la  a  u^-cycle  of  tha 

0  0 
procaaa  v  ,  aa  vail  aa  tha  sum  of  cha  paaaaga  ciaaa  la  a  u  -cycle,  dcas 

not  dapaad  on  cha  Idancldas  of  cha  jobs  la  successive  configurations  of 

tha  job  stack.  It  follows  chat  va  can  work  vich  tha  procaaa  W«\V(t) ; tiOx 

daflaad  by 

WCt>  -  (K(t),2(t))  . 

Here  He'1  ccrraaponds  to  tha  liaaar  job  stack,  and  Kvt)  is  tha  last  state 
visited  by  tha  Markov  chain  £«{Z(t):tiO}  bafora  Jumping  to  2  vO  •  Tha 
procaaa  V  is  an  irraduclbla,  positivt  racurranc  continuous  tiaa  Markov 
chain  vich  a  stata  spaca  that  is  a  subsat  of  2*2.  Moca  that  in  gar.aral 
tha  stata  spaca  of  V  is  much  smallar  than  that  of  V^,  and  that  working 
with  cha  procass  V  is  computationally  advancagacus . 


Tha  computations  rase  on  cha  dafiaition  of  two  particular  functions 
(f  and  g)  daflaad  on  tha  stata  spaca  of  W  and  taking  valuas  in  tha 
sat  {0,1}.  Va  dafina  tha  functions  f  and  g  in  tans  of  functions  f  ^  and  gv 
dafinad  on  F “ ,  tha  stata  spaca  of  tha  procass  V1*.  Va  cake  tha  function  f'* 


Va  taka 

to  ba  tha  indicator  function,  1  of  cha  sac  j1"  which  defines  tha  starts 


J 


ot  ;»M{«  tiaes  irrespective  of  job  identity,  i.«.f  ;ar 
(r.a^, . . . ,a^,a' ,a^ . a^)«F°, 

U'°- . V!''3i . =*>  ’  :sot!’“i . . «»>  ■ 

rSu<  ~  4  passage  diae  far  some  job  starts  vh«n  V°  hits 

U,31 . V2’’ai. . 9*  :h#n  f°-:*  *Vc=*  for  each  (*,*•)  ia  :ht 

'  ^•5w  ai»  •  •  •  »a}i.a1.  •  •  • .  and  n^  such  that 

(*,al . V,,*ai . 9<?0*  ?or  *uch  *  C*,z'),  ve  define 

f(X’Z,)  "  f0(Z,al . V*’’ai . 9  •  (5.5 


The  ruaccian  i  is  vail  defined  since,  for  fixed  z 
is  independent  of  its  other  arguments. 


ana  z  ,  the  function 


.0 


For  an  element  (c,n. 


•  9«?0,  :: 


he  value  of  the 


function  i J  is  tne  number  of  passage  tines  that  start  or  are  under--, 
the  configuration  of  the  job  stack  is  o’.  Formally,  for 
v*,a1,...,aJ(1s,,a^ . 9«f°,  ve  define 

*  ‘2,al . V2’  ,al . 9  *  £  i-Chu\k))  . 

k»l  * 

hen,  for  ,c,c’)  ia  the  state  space  of  V,  ve  set 

’  **<*'*! . V*"9"-'9  • 


va  •  0/ 


•aaSii: t4«i«  ^ir.g  the  process  V  is  that  the  number  of  cussaoe 


•3^3  ^a*a'i  V  «*  ~  £4  - 


23 


V1 


£  f<V  . 

j-0  J 


*“  ""  “  *  “ . —  -  -  ot  2o  u 


r 


«(W(*))di 


'*i  [r«ap.  k^]  ia  tht  tin#  of  ch^  flrjt  r#t  f 

Ch*  JUOP  °h4in  <V*»I  to  th.  fixed  scat.  w*  ’  '  '  Pr°C*SS 

w  .  xhe  r.curn  .pjf. 


Chen  v°.(X(lV 


0 

a  •  (  9  n  0 

. . V2  *ni . “J> 


By  direcc  application  of  Ei 


<lu* Cions  0.1)  end  (3.2),  v.  h4V. 


PROPOSITION  5.1.  Lac  f  b. 
cha  cransicio 


<“nctloa  *  -^«;on  a. 5),  a 

-““U  *'  " *rr"  r,„ 


Et'Y.(f)}  .  E 


and 


(k-0  *  ) 


1  ) 

((W  -  u-.,a)_ie 


E{(Vf))2}  -  (i-QR)-ih  t 

hVh*r*  <l  **  th*  ““  of  th-  “*«  —  =0  ch.  „o  ta 


and 


J 


Jt  . 


:6 


? tom  Propoaicioo  (5.1)  va  obcain  I  j(S^:  and  £  ,{(*?)*}  according  co 

u  uu 

£  Q(JCJ)  -  £  0(Mf)} 


(5.7) 


And 


£  •  £  .{(Y  (i))2}  . 

w 


(3.3) 


3y  diracc  appiicacion  oi  Iquaciona  v3.3)  And  (3.1),  va  Sava 


PROPOSITION  5.2.  Lac  j  ba  cha  iunccion  iaiinad  by  Tquacion  ;3.a),  And  R 
ba  cha  Jump  aaerijc  And  q  cha  '.’actor  o£  raca  oaraoacars  ior  holding  ciaaa 
In  cha  concinuoua  ciaa  Mar’acv  chain  Than 


\  r;’  > 

!lV*,;  •  S)J,  ' 


vl- vR)”^<.g»q~*)  , 


and 


SUVi))*:  -  ,I-0.V  “ \>hjq 


-1  -L 


whara  ca  cha  ciaa  oi  cha  firsc  racurn  co  cha  jcaca  vj  in  V,  and 

Sa ’ * a  ?  .  V  • \  \ 


Va  uaa  chi j  rasulc  co  obcain 


and 


.5.9) 


.5.10’' 


I 
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Using  che  r«do  formula,  Equations  (5.7)  and  (5.9)  yield  cha  valua  of 
r^.  Analogous  to  Proposition  4.3  we  hava 


PROPOSITION  5.3.  For  EfY^f)}  and  EtY^g)}  given  by  Propositiona  5.1  and 
3.2, 


l)  I  1  *(W(s))ds  £  *<V 

ro  j-o 


) 

r 

where  h-uVVEtY^OW^Y^gjMg'q’Vf  • 


(I 


Wa  usa  this  rasult  to  obtain 


and  an  axprassion  for  Che  varlanca  constant  follows  from  Lquacions  (5.3), 
(5.10),  and  (3.11). 


6.  NUMERICAL  RESULTS 

Wa  considar  tha  closad  network  of  quauas  of  Figure  1  and  the  limiting 
passage  times  P  and  R  therein.  Upon  completion  of  service  to  a  Job  at 
center  1,  in  accordance  with  a  binary-valued  variable  y,  tha  Job  Joins  the 
tail  of  cha  queue  la  center  1  (when  y*l)  or  Joins  the  tall  of  che  queue  in 
center  2  (when  y«0) .  Neither  center  1  nor  center  2  service  is  subject  to 
interruption.  We  assume  that  both  centers  provide  service  according  to 
a  FCFS  (first-come,  first-served)  queueing  discipline.  The  limiting 
passage  time  ?  is  associated  with  the  time  measured  from  entrance  by  a  Job 
into  che  center  1  queue  upon  completion  of  a  center  2  service  until  the  Job 
next  enters  the  center  2  queue.  Similarly,  the  response  time  R  is  associated 


:s 


wish  tha  tiaa  bacvaaa  succaasiva  aatrancas  of  *  Job  into  tha  cantar  1  quaua 
upon  comp  laden  of  *  cantar  2  sarvlca. 

Va  maka  tha  assumptions  that  (i)  sarvlca  tiaas  at  cantar*  1  and  2  art 
mutually  ladapandant  and  Cii)  succaasiva  sarvlca  daas  at  cantar  1  fora  a 
saquanca  of  1.1.4.  random  variablas  axpcnantially  distributad  with  rata 
paramacar  \  ,  1*1,2.  In  addition,  va  as sum#  that  tha  routing  varlabla  v  la 
a  Barnoulli  random  varlabla,  and  valuta  of  J  at  succaasiva  oomplatlona  of 
sarvlca  at  oantar  1  fora  a  saquanca  of  l.l.d.  random  variablas,  ladapandant 
of  tha  sarvlca  tiaas  at  cantars  1  and  2.  Va  dances  by  p  tha  probability 
that  tha  routing  random  varlabla  v  takas  tha  valua  1. 
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Similarly,  the  sets  and  defining  the  terminations  of  the  peeaege 
time  P  ere 

Bl  -  { (1,1) :0<1sn)  , 


end 


B2  -  { (1-1,1) :0<1SN}  . 

For  the  response  tlae  R,  the  sets  A^  end  A,,  ere  the  sene  es  for  the  pessege 
tine  P,  but  B^-A^  end  B^-A^. 


In  connection  with  the  merited  job  method,  the  continuous  tine 
process  V»{ (L(t) ,X(t)) : tkO),  where  L(t)  Is  the  lest  state  visited  by  the 
Markov  chain  JC  before  jumping  to  X(t) ,  has  state  space  F,  where 

F  -  {(i,j,i+l,j+l):Osi<N,  IS J<N}  u  { (i,N, 1+1,1) :0si<N}  u 

{ (l.J  »i“l» j ) * 0<1SN,  lSjSN)  u  {(1,1 ,1,1) :1<1SN}  . 

The  subsets  of  F  defining  the  starts  and  terminations  of  passage  times 
for  the  marked  Job  are 


S  -  {(i,N, 1+1,1) :0si<N} 

and 

T  -  {  (1,1, 1-1,1)  :0<isN}  . 

Tables  1  and  2  give  theoretical  values  for  simulation  of  the  closed 


network  of  queues  by  the  marked  Job  method.  Numerical  results  are 
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diaplayad  for  eh*  naan  of  eh*  raapoaa*  ela*  %  and  corr**pondiag  raaulca 
for  eh*  pasaag*  eia*  ?  tri  ia  partach**** .  For  eh*  cat*  of  S»2  Job* 

(Tabi*  1) ,  eh*  «*e  S-{ (0,2,1, 1) ,  (1,2,1, 1)>.  Wich  \^L,  \2-0 .5 ,  and  p-0.75, 
eh*  auaarlcal  rasules  show  chac  oa  eh*  average  0-cycl*a  d«fia«d  by  r«cura* 
eo  eh*  scat*  (0,2,1, 1)  *r*  rwic*  ««  long  ««  ehoa*  d«fla«d  by  eh*  r« curas 
eo  eh*  seae*  (1,2, 1,1).  Moe*  ehae  a*  axp«ce*d,  eh*  quantici**  a 
(a*  wall  a*  (SgCo^f)^ ar*  eh*  saa*  for  eh*  eve  r*cura  seae**. 
Tabl*  2  give*  results  for  job*.  H*r*  char*  ar*  four  pcaalbl*  raeura 
seaeas,  and  for  eh*  parameter  valuaa  s*l«ce«d,  return*  eo  eh* 
jeae*  (3, 4,4,1)  occur  aoae  frequently,  and  oa  eh*  *v«rag«  «ighc  cia*s  nor* 
ofc*a  ehaa  raeura*  eo  eha  «eae*  (0, 4,1,1). 


V*  now  turn  eo  eh*  decomposition  aaehod.  Th*  process 


.,0 


-  { (2(e)  ,S^(e) ....  ,if*(e)) : eaO} 


ha*  seae*  spac*  £  ,  vh«r« 


>0  _  , 


2  •  v  (i.a^,  . . .  ,a^)  :0siSN;  Isa^ . a^«J ;  a^tt,  for  i*J  ; 

7h*  underlying  conciauous  eia*  process  :  eaO:  d*£ia*d  by 


vC(c)  •  v*-0vc),:c°(e))  , 


vh«r*  L1' (e)  1*  eh*  la«e  seae* 
eo  XJCe),  ha*  seae*  spac*  F°. 


'/is lead  by  eh*  Marfcov  ehaia  ^  b*for*  Jumping 
Th«  subsae*  of  Jw  defining  seares  and 


in* ties*  of  passage  eio*s  ar* 
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T  ■  ( ( 1  , , rLy •  •  •  t^|)  •  <X tsN  { 

lsaj«njsN,  for  lsjsN;  o^dn^  for  kdj}  . 

Tha  procasa  J«{2(t) :tiO}  has  stitt  spaca  D«{0,1, . . . ,N},  cha  sac  D°*(0}t 
and  cha  sac  1'^  daflnlng  cyclaa  of  cha  procasa  V®  is 

U°  ■  {(l,n^, . . .  .o^.O.n^, . . .  ,a^) :  lsn^, . . .  ,rVjjSM;  n^dn  for  l*j}  , 

Tha  scaca  spaca  of  cha  acochaacic  procasa  V  la 


{ (1,1+1) :0ilsN-l}  u  {(1,1-1) :1S1SN}  , 
and  cha  scaca  w^«(i,0). 


Tabla  3  givas  thaotacical  valuaa  for  simulation  of  tha  olcsad  aacverk 
of  quauas  by  cha  dacomposicion  aachod  for  cha  quantity  £;?}.  Tha  tabla 
glvaa  rasulta  for  N*1  to  N‘«<*  Jobs,  and  tha  paramatar  valuaa  art  cha  sama 
aa  In  Tablaa  1  and  2.  For  h’«2  Jobs,  cha  valua  of  cha  quantity  a'\u^)»av'1 
which  aaasuraa  cha  statistical  afflclancy  of  tha  dacompoaltion  machcd  la 
16.5-6.  Tha  corraspondlng  valua  from  tabla  1  for  tha  aarkad  Job  aachod  is 
20.390.  thus,  for  chaaa  paraaacar  valuaa  tha  dacompoaltion  aathod  Is 
approxiaacaly  21  parcanc  acra  afficiant  than  tha  aarkad  Job  aathod.  For 
N*4  Jobs,  tha  dacoapoaltion  aathod  la  -1  parcant  aora  afficiant. 


i 


J 
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Muaerical  results  pertaining  Co  eh«  statistical  efficiency  of  the 
decomposition  as c hod  in  simulation  of  eh*  closed  necvork  of  queues  appear 
in  Table  4.  Tor  N«1  to  M«b  jobs,  th*  cable  give*  theoretical  value*  of 
Che  quantities  r°  and  e°  for  three  secs  of  parameters  values.  W*  hold  the 
value  of  \j»l  and  p-O.TJ  fixed,  but  vary  Table  5  give*  a  comparison 

of  ch*  statistical  efficiency  of  the  narked  job  and  decomposition  aethcds 
for  the  same  sets  of  parameter  values. 
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TABLE  1 


ThaoracicAl  Valuta  for  Ch«  Harkad  Job  Machod. 
P«a a Aft  TIba  R  CP)  La  Cloaad  Nacvork  of  Quauaa. 
N"2,  Xj«l,0,  Xj-0.5,  p-0.75. 


PArABACtr 


W 


I  Mi  ) 


W 


0  /E..VM, } 
0  1 


Raturo  Scat#  of  V«{V(t):ciO} 


\  1  /  2 

VP  ff/t0{Ml> 


— 1 

CO, 2, 1,1) 

<1.2, 1,1) 

24.0 

12.0 

23.0 

14.0 

C20.0) 

(10.0) 

3.0 

1.5 

C3.0) 

(1.5) 

9.333 

9.333 

v.6.667) 

(6.667) 

140.26: 

70.133 

(129.067) 

^6-. 533) 

>  % V 

46.756 

*6.756 

C43.022) 

(43.022) 

20.390 

20.390 

(20.033) 

(20.033) 

TAALS  2 

Theoretical  V«lu««  for  the  Marked  Job  Mac hod . 
?uit|«  Tiae  ?  (R)  ia  Closed  Rework  of  Queues. 
M«4,  X^«1.G,  Xj-0.5,  ?*0.75. 
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TABLE  3 


Thaoratical  Valua#  for  tha  Dacomposic: 
Paaaaga  Tia«  P  la  Cloaad  Nacvork  ol 
.  Aj«1.0,  X^-O.S,  p-0.73. 


Parana car 

N-l 

N-2 

1  . . .  . 

N-3  N-4 

E  o{a?} 

6  a  G 

24.0 

30.0 

62.0 

u 

i  K?  ol 

*„Z  p? 

4.0 

20.0 

63.0 

196.0 

u  (j-1  J) 

i 

E  0{kJ} 

1.0 

3.0 

7.0 

15.0 

u 

K1  n)  - 

E0  £ 

4.0 

6.667 

9.714 

13.067  | 

u  ( j»L  J  )  u 

1 

»:>• 

16.0 

176.0 

1023.673 

i 

1 

4317.227 

<V Z/s  otKi; 

16.0 

53.667 

146.249 

237.315 

u 

l 

(E  o(ar)l/2  4'e  o"?> 

9.798 

16.546 

25.035 

34.491 

\  u  /  u 

1 

_  ___ 

_  _ 

36 


TABLE  4 

Statistic*!  Efficiency  of  the  Decomposition  Method. 
Passage  Tiae  ?  in  Closed  Network  of  Queues. 


?  ■ 

— — — — i 

0.75 

■9 

?  • 

0.75 

\  • 

1.0 

K 

Vl  * 

1.0 

0.125 

\i  - 

0.5 

* 

*• 

0 

0 

0 

— 

0 

o 

0 

r 

e 

r 

* 

r 

e 

^  %  0 

13.356 

4.0 

12.31* 

-.0 

9.793 

3.333 

19.956 

6.0 

17.664 

6.667 

16 . 5-6 

6.236 

27.330 

3.0 

26.123 

9.714 

25.035 

6.933 

35.139 

:  o.o 

36.606 

13.06* 

3-. -91  , 

7.335 

42.597 

12.0 

-9.107 

16.6-5 

—  .296 

7.619 

•*9.063 

,  , 

63 . 6-5 

20.331 

5-. 021 
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APPENDIX 

This  appendix  is  davoetd  to  a  proof  of  -ha  result  given  in 
Equation  (3.5).  Recall  chat  X-{X(t):tiC}  is  a  continuous  state  Markov 
chain  vich  finite  state  space  E*{0,1,2, . . . ,N; .  We  assume  £  is  irreducible 
and  hence  positive  recurrent.  The  mean  holding  time  in  state  i  is  qv *  and 
the  embedded  jump  chain  {X  :nifl}  has  transition  matrix  R.  Again  f  and  g 
are  real-valued  functions  defined  on  E  and  we  view  them  as  column  vectors; 

e.g.,  f»(f(0) . f (N) ) .  Ve  take  0  as  the  return  state  and  as  before  let 

3,  denote  the  time  of  first  entrance  of  1C  to  3  and  a.  denote  the  time  of 

1  -  i 

first  entrance  of  (X  :niO;  to  3.  lac  X(0)»C,  7  «0,  and  7  (nil)  be  the 

3  U  a 

random  time  at  which  the  r.ch  transition  of  X  takas  place.  Since  X  is 
irreducible  and  E  is  finite,  T  ••  with  probability  one. 


PROPOSITION.  Under  the  conditions  aoeve 


i  c 3.  v1  ) 

E  *  r\X(s))ds  £  5(X  )  -  C-,R)  ~n  , 

I  *'n  --0  3  1 


where 


h  ■  (I*cR)”i,3>"f >q  “  V  i(Ir-R)  “(f»q  “)  j  »g.  -  f»o  *»* 
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Proof.  First  w«  dacomposa  cha  axprassion  in  Equation  (A.l)  as  follows: 

ri  *  V1  V1  V1 

1,  t(x(,)>d'&  •“»>■£  "vvvt  .<v 

'  n?0  lo  f(X»Urn»l‘T»,1{«1>n)‘<X-)1(5l>») 

■  E  E  ■„».  ♦  E  E*b  -  f  .  b  ,  (A.; 

n-0  m-0  n  a  n-0  m-n  n  m  n-0  an 

whsrs  sq  and  b^  art  dafinad  la  cha  obvious  way.  Sow  taka  axpac cations  on 
both  sidas  of  Equation  (A. 2).  Wa  computa  tha  chraa  axpactations  on  cha 
right-hand  sida  of  Equation  (A. 2)  saparataly.  Tha  incarchangas  balow  of 
Ej_  and  Z  can  ba  Juscifiad  sinca  cha  procass  JJ  is  poslciva  racurrant  and  tha 
stata  spaca  is  finica.  First  tna  assy  cans: 

j  •  )  •  .  V 

E,  Z  Oal-  Z  E.  fCX  >  1<T  .-T  )( 

M  n-0  n  n  '  ns<2  M  n  n  1  o  ^>n;  n*l  n  J 


lH&  •‘vj  • 


From  Equation  (3.1)  wa  can  concluda  chat 


_ J 
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*!S  v.l  •  <i-aw:l«*rli««)  . 


(A. 3) 


Th«  first  tara  oa  tha  right-hand  aid*  of  Equation  (A. 2)  ta  handlad  In  a 
similar  fashion.  Aftar  conditioning  on  {XjiJSn}  as  abova,  ua  hava 

EmIo  £  *“S*1  ‘  £  £  Et!f0t»),"laaUli1>i.:'‘a»U{S1>«>! 

•  E  t  vW  • 

n-0  a-0  m 

Mov  incarchanga  tha  sums  and  aaka  a  chang a  of  variablaa  to  obtain 

h\  E  t  s»,j  -  E  E  lt(eaV 

( n-0  n-0  )  a-0  o-a 


Z  Z  V;Wa} 

a-0  n-0 


Z  Z  Z  Z  &  0r?ki(J)JOOq~1\iO  . 

a-0  n-0  j«E  k<E  u  J  ^ 


Siaca  T  ,Ra-(I-.arL,  wa  concluda  that 
n~0  J  0 


Tha  aacond  tarts  on  tha  right-hand  aids  of  Equation  y.A.2)  raquiras  i  scmawna 
dlffartat  tonditionisg  arguaant. 
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characteristics  of  "passage  times"  in  certain  closed  networks  of  queues. 
These  estimation  procedures  are  based  on  the  regenerative  method  for 
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traverse  a  portion  of  the  network.  Such  quantities  are  Important  In 
computer  and  communication  system  models,  and  In  this  context,  quantities 
other  chan  mean  values  are  of  interest.  From  a  single  simulation  run, 
our  passage  time  simulation  methods  provide  both  point  estimates  and 
confidence  intervals. 

We  consider  here  the  calculation  of  variance  constants  entering  into 
central  limit  theorems  used  to  obtain  confidence  intervals  from  passage 
time  simulations.  Using  results  of  Hordljk,  Iglehart,  and  Schassberger 
for  Che  calculation  of  moments  in  discrete  time  and  continuous  time  Markov 
chains,  we  calculate  variance  constants  pertinent  to  mean  passage  times. 

We  do  this  first  for  the  "marked  job  method"  for  passage  time  simulation 
which  is  based  on  Che  tracking  of  a  distinguished  job,  and  then  for  the 
"decomposition  method"  in  which  observed  passage  times  for  all  of  the  jobs 
enter  into  the  construction  of  point  and  interval  estimates.  The  results 
of  this  paper  provide  a  means  of  comparing  the  statistical  efficiency  of 
the  two  estimation  methods. 
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